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the competing order enhance the phase-space available for low-energy quantum fluctuations and so 
self-consistently lower the free energy. We demonstrate that this quantum order-by-disorder mech- 
anism leads to instabilities towards the formation of spiral and d-wave spin nematic phases close to 
itinerant ferromagnetic quantum critical points in three spatial dimensions. 
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I. INTRODUCTION 

A variety of unusual phenomena that do not fit into 
the framework of Fermi-liquid theory have been observed 
in the vicinity of quantum phase transitions. An inter- 
esting example is the emergence of new phases near to 
putative quantum critical points, observed in many ex- 
periments. Examples include the onset of superconduc- 
tivity close to the i tine rant ferromagnetic quantum crit- 
ical point of UGe 2 f^l an anomalous anisotropic phase 
around th e m etamagnetic quantum critical end point of 
SrsR^OvpIAIa possible inhomogeneous magnetic state in 
ZrZn 2 por the unusual partially ordered phase of MnSi. 6 
This has led to the speculation that the onset of new 
phases might represent a generic principle.^ 

In his pioneering work, HertzSl studied the paramagnet- 
to-ferromagnet quantum phase transition of itinerant 
fermions that occurs by varying the exchange coupling 
between electron spins. He derived an effective action for 
dynamical fluctuations of the bosonic order parameter. 
Later, MilliJ^ used this approach to calculate tempera- 
ture dependencies of the correlation length, susceptibility 
and specific heat. In the past decade several authors^HIl 
carried out diagrammatic calculations that extend be- 
yond Hertz-Millis theory. They showed that the free en- 
ergy contains a non-analytic dependence on the order 
parameter and its gradients. It was arguecP^ that non- 
analytic terms occur due to additional soft particle-hole 
modes that couple to the order parameter. These non- 
analytic corrections can render the transitions weakly 
first order at low temperatures and lead to the insta- 
bility of quantum critical points to the formation of new 
phases. 

We present an alternative, analytical approach - quan- 
tum order-by-disorder - which is able to predict new 
phases and provides an intuitive physical picture of the 
problem. Moreover it results in relatively simple calcu- 



lations, accessible to a broad audience. It relies on the 
idea that certain deformations of the Fermi surface en- 
hance the phase space available for quantum fluctuations 
and thus self-consistently lower the free energy. This re- 
sults in new phases near to the putative quantum crit- 
ical point. Note that in the familiar realizations}^^ 
the quantum order-by-disorder mechanism is driven by 
bosonic order-parameter fluctuations. Here, the under- 
lying fermionic statistics of the excitations and Pauli 
blocking of phase space is very important. Recently, it 
has been demonstrated^ that the formation of the enig- 
matic partially ordered phase of MnSP with its pecu- 
liar magnetic-ordering pattern can be explained by the 
fermionic quantum order-by-disorder mechanism. 

Our work builds upon the work of Conduit et a/p2l in 
which a spiral phase was predicted close to the itinerant 
ferromagnetic quantum critical point in three spatial di- 
mensions. In that work, the quantum order-by-disorder 
approach was used with a numerical evaluation of the 
fluctuation corrections to the free energy in the presence 
of a spiral state. We develop an analytical approach that 
ultimately allows us to extend the framework to include 
new phases such as the spin nematic, a Pomeranchuk- 
type instability in which the net magnetization is absent 
and the Fermi-surface deformations for spin-up and spin- 
down electrons are of opposite signpH23 

This is achieved through a Ginzburg-Landau expansion 
in the vicinity of the finite temperature tricritical point. 
We calculate closed expression for the Ginzburg-Landau 
coefficients of a uniform ferromagnet and evaluate them 
analytically in the limit T — > 0. The Ginzburg-Landau 
coefficients of the spiral ferromagnet are related to the 
coefficients of uniform ferromagnet by averages of certain 
angular functions. 

For the spin nematic we first develop an expansion of 
the generating function in powers of the field conjugate to 
the spin-nematic order parameter. Similarly, the coeffi- 
cients in this expansion are related to uniform Ginzburg- 
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Landau coefficients by averages of certain angular func- 
tions. A Legendre transform of the generating function 
recovers the expression for the free energy. By including 
small deviations from the isotropic free-electron disper- 
sion, we are able to obtain phase diagrams relevant to a 
broader range of experimental systems. 

This paper is organized as follows: In Section II we 
discuss the key ideas of quantum order-by-disorder and 
outline its mathematical setting. In Section III we pro- 
ceed to construct the Ginzburg-Landau expansions for 
the uniform ferromagnetic, spin-spiral, and spin-nematic 
states. This enables us to construct a phase diagram in 
Section IV. Finally, in Section V we summarize our re- 
sults and suggest directions for future work. 



II. QUANTUM ORDER-BY-DISORDER 

The central idea of quantum order-by-disorder is that 
certain phases have more low-energy quantum fluctua- 
tions associated with them. This lowers their free energy 
and renders them stable. The effect of the lowering of 
the free energy already becomes evident in second order 
perturbation theory; the second order contribution to the 
free energy of the ground-state is always negative. The 
mechanism is similar in some ways to the entropic low- 
ering of the free energy in classical systems. Equivalent 
results for the contribution of quantum fluctuations to 
the free energy can be derived starting from a functional 
integral approach. 

A well known example of quantum order-by-disorder is 
that of a quantum antiferromagnet. If the electron spins 
are oriented ferromagnetically, no virtual electron hop- 
ping is allowed due to the Pauli exclusion principle. On 
the other hand, when the spins are antiparallel, electron 
hopping is allowed. This hopping lowers the free energy 
of the system through second order perturbation theory. 
In this way the antiferromagnetic phase is stabilized due 
to the effect of quantum fluctuations. 



A. Perturbation Theory 

Let us now begin to develop this general heuristic pic- 
ture into an explicit calculation. Our starting points is 
the free electron system in d = 3 spatial dimensions in- 
teracting through Hubbard point repulsion 
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Here — ^- is the isotropic free-electron dispersion, 
fi denotes the chemical potential and n±(r) densities of 
spin up/down electrons. Note that later on, we will in- 
clude small anisotropic deformations to make the disper- 
sion more tight-binding like. The strength of the contact 
interaction is given by g. The mean-field free energy is 



given by 

^mf = -i^ln(l + e-^-^)+ 5 J d 3 rAf 2 (r)(2) 

k,er 

where j3 represents the inverse temperature, e£ is the 
mean-field dispersion in the presence of certain type of 
order, and M(r) is the magnetization vector. In this work 
we will not consider phases with spatial charge modula- 
tions. 

The effects of fluctuations are already seen in self- 
consistent second order perturbation theory. The fluc- 
tuation corrections to the free energy are given by 

n n t U Z (1 — n t )(1 — nZ ) 

j- fl = -2 5 2 + i - _ + _ - — . ( 3 ) 

k 4 ...k 4 e k!+ e k 2 6 k 3 £ k 4 

where the prime over the sum indicates momentum con- 
servation, + k 2 = k 3 + k4, and for brevity, we have 
written the Fermi functions as 

<:=n( £ £) = (e^-^ + l)- 1 . (4) 

Note that the fluctuation corrections to the free en- 
ergy are calculated self-consistently; the energies enter- 
ing the Fermi functions are the mean-field dispersions in 
the presence of a given type of order. 

From Eq.([3| we see that the fluctuations correspond 
to excitations of virtual pairs of particle-hole pairs of op- 
posite spin and equal and opposite momenta (spin up 
particle-hole pairs carry momentum ki — k 3 and spin 
down particle-hole pairs carry momentum k 2 — k4). Since 
we need to put in energy to create the particle-hole pairs, 
the denominator of ([3| is always positive, which results 
in negative contributions to the free energy. Certain de- 
formations of the Fermi-surface enhance the phase space 
available for low energy, virtual particle-hole excitations 
and in that way self-consistently stabilize new phases. 
Ferromagnetic, spiral or spin-nematic Fermi-surface dis- 
tortions which are shown schematically in Fig. (JlJ all 
open up extra phase-space for the low-energy particle- 
hole pairs to form. 

Careful inspection of Eq. ^ reveals that the term 

3r = -2 5 2 E + - n ^ k ; _ _ (5) 

k 4 ...k 4 +£ k 2 C k 3 £ k 4 

contained in gives an unphysical divergent contribu- 
tion to the free energy. To avoid this, we need to take 
into account the renormalization of the interaction ma- 
trix element g^ We perform a self-consistent perturba- 
tive expansion around a mean-field solution. Instead of 
using the momentum to label the eigenstates of the free- 
electron Hamiltonian, we use it to label first-order shifted 
states, 

|kt,U) = |kt,U)o 

V o(Pt,ql|H i " t |kt i U) 

+ 2^ + - _ + _ - — IP 1)3-1-/0 

p ^k,q#i e k + e i £ p £ q 

(6) 
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FIG. 1: (Color online) Distortions of the Fermi surface 
(dashed lines) enhance the phase-space for quantum fluctu- 
ations, (a) uniform ferromagnet (b) spiral and (c) d-wave 
spin nematic. Quantum fluctuations correspond to excita- 
tions of pairs of particle-hole pairs of opposite spin and equal 
and opposite momenta. 



where |k f, 1 l) labels the two-particle free electron 
state, |k |, 1 I) labels the first order corrected two- 
particle state, e£ are taken self-consistently to be the 
mean-field electron dispersions and W nt represents the 
interaction Hamiltonian. With this identification, we 
must also make a corresponding alteration to the inter- 
action strength, 



k 3 ,k 4 k i 



(7) 



This renormalization corresponds to a one-loop cor- 
rection to the interaction strength and leads to a regular 
expression for the free energy, 



'+n k2 (n : 



ki...k 4 



(8) 



B. Functional Integral Derivation 

We next sketch how the sa me res ult can be derived 
using field theoretical methods! 20 * 25 ! This approach re- 
veals immediate connection between the second order 
self-consistent perturbation theory, outlined above, and 
field theoretical calculations that explicitly show non- 
analytic behavior of the free energy. We start from the 



fermionic partition function, 

Z = J V($,i))e- S &^\ 

S$, V>] = jf " dr J d 3 r + H$, i>)] , (9) 

where ip = (i[> + ,ip_) T an d V> = (i/>+,V>_) denote Grass- 
man fields which vary throughout space and imaginary 
time, and the Hamiltonian is given by Eq. Q. After per- 
forming a Hubbard-Stratonovich decoupling of the inter- 
action in spin (<p) and charge (p) channels we obtain 

Z = J V{^,%P)Vct)Vpe' s ^^^P\ 

S[i>^,<t>,p] = J ^(Gu 1 + g(p - 4> ■ a))i> 

+9 f(d> 2 -P 2 ), (10) 



where Gq 1 is the free-electron Green function and er de- 
notes the vector of Pauli matrices. Integrating out the 
fermions, we obtain 

Z = J V(j}Vpe- s[<l, ' p \ (11) 
S[<M - -Tr In [<V + <?(/>-<'- -0)1 +9 \{<?-p 2 )- 



So far, all the steps are the same as in Hertz-Millis 
theory. However, in that case the aim was to derive an 
effective action for dynamical fluctuations of the bosonic 
order parameter in the paramagnetic state. In contrast, 
we wish to derive a Ginzburg-Landau expansion in the 
static order parameter. In order to do this we separate 
4> and p into zero-frequency (po, M) and finite-frequency 
parts (p,(/>): p = po + p and cjy = M + 4>. The action then 
becomes 



SI 



-Tr In 



G Q 1 + gpo - gtr.M + g(p - <j> • tr) 



+ g J (M 2 + ^-p 2 ). 



(12) 



We expand this action to quadratic order in finite- 
frequency fluctuations and integrate them out, yielding 
the following expression for the free energy: 



T[M] = gM 2 - TrlnCT 1 
+ -Trln(l + 2gYY ] 



ff 2 n+-n- 



1 



Trm(i-^n ++ n--) . 



(13) 



■7" a, i 



The term Fg. ± arises from transverse magnetic fluctu- 
ations, while the term Jj n accounts for charge-density 
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and longitudinal magnetic fluctuations. The polarization 
bubbles are given by 

IT"'(q,w) = ^^^(k-qV-ij), (14) 

k,a/ 

where G" 1 = G^ 1 +gpo—agM . The polarization bubbles 
W 7 explicitly depend on the interaction strength g. 

Next, we perform a rather strange expansion in g. 
We expand the Trln-terms to second order in g by only 
expanding in powers of g that stand in front of the 
polarization bubbles, while keeping the full g dependence 
of the polarization bubbles as it is. This looks like a 
second order expansion in g, but self-consistency actually 
implies resummation of certain classes of contributions 
to infinite order. This expansion captures the relevant 
physics, as we will see later on. After performing the 
summations over Matsubara frequencies, we arrive at 
the expression Further, we need to renormalizc g, to 
cancel the ultraviolet divergence. Doing so according to 
Eq. (JtJ) recovers expression Eq. (|Sj) for the free energy. 

In summary, quantum order-by-disorder provides a 
physical picture for the formation of new phases due to 
quantum fluctuations. As we will see later, certain de- 
formations of the Fermi surface enhance the phase space 
available for quantum fluctuations and in that way lower 
the free energy. This is already apparent in second order 
perturbation theory and can also be derived from a func- 
tional integral approach. Next, we want to expand the 
free energy in powers of the order parameter (which en- 
ters through the mean-field dispersion) . This will enable 
us to construct the Ginzburg-Landau expansion and to 
analyze the phase diagram. 



III. GINZBURG-LANDAU EXPANSION 

We wish to determine the phase diagram of the near 
critical itinerant ferromagnet allowing for the generation 
of new phases near the quantum critical point. In order 
to obtain the phase diagram, we develop a Ginzburg- 
Landau expansion of the free energy in powers of the 
order parameters for the various types of phases that 
might form. The expansion is valid in the vicinity of the 
tricritical point (explained later), where the value of the 
order parameters is sufficiently small. 

It turns out that our task is simplified considerably by 
relationships between the expansion coefficients for the 
different types of order and those for the uniform ferro- 
magnet. We begin therefore with an explicit evaluation 
of the coefficients for the uniform ferromagnet. 

Next, we allow for spatial modulations of the ferromag- 
netic order; in particular we consider a spiral modulation 
of the magnetization. We use the fact that the free en- 
ergy can be expressed (to all orders) as a functional of 
the mean-field electron dispersion in the presence of the 



spiral order. We show how the coefficients of the ex- 
pansion in the spiral ordering wave vector Q are related 
(by angular averages of certain functions) to those of the 
uniform ferromagnet. 

For other order parameters, not driven in mean-field, 
we introduce a field conjugate to the order parameter 
and construct an expansion of the generating function in 
terms of the conjugate field. We are able to relate the 
coefficients of the generating function to the Ginzburg- 
Landau coefficients of a uniform ferromagnet by per- 
forming some simple angular integrals. We use a Leg- 
endre transformation of the generating function to re- 
cover the Ginzburg-Landau expansion. Quantum fluc- 
tuations generate an interaction in the new channel. A 
similar mechanism is well known in spin-fluctuation the- 
ory where s uperc onductivity is mediated through spin- 
fluctuations I2SED j n w hat follows, we will concentrate on 
the case of a d-wave spin nematic. 

Finally, we allow for a more generic energy disper- 
sion by considering small anisotropic deviations from the 
isotropic free-electron dispersion. We calculate the cor- 
rections to the coefficients of the Ginzburg-Landau ex- 
pansion due to the anisotropic distortion. The coeffi- 
cients of this expansion are proportional to parts of the 
Ginzburg-Landau coefficients of the uniform ferromagnet 
in the presence of an isotropic dispersion. The propor- 
tionality factors are determined by angular averages of 
functions that encode the specific form of the deviation 
from the isotropic free-electron dispersion. 

A. Uniform Ferromagnet 

The dispersion of the uniform ferromagnet is given by 
e CT (k) = k — agM. We Taylor expand the free energy in 
powers of M, 

F[M) = aM 2 + f3M 4 + 7 Af 6 + (15) 

where the Ginzburg-Landau coefficients are functions of 
interaction strength and temperature, a — a(g,T), and 
similarly for /3 and 7. 

1. Mean-Field Coefficients 

The expansion of the mean-field free energy Eq. ^ in 
powers of M leads to the following coefficients: 

umf = .g + .g 2 ^n (1) (e k ), 

k 

0MF = !.9 4 £" (3 W 
k 

IMF = |S 6 5> (5 W (16) 
k 

The remaining integrals over derivatives of Fermi func- 
tions are straightforward to compute for the present k 2 
dispersion. 
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2. Fluctuation Contributions to Coefficients 

Here we outline the main steps in the calculation of 
the fluctuation contribution to the Ginzburg-Landau co- 
efficients. The detailed calculation is given in the ap- 
pendix. We start by writing the fluctuation corrections 
to the free energy in terms of modified particle-hole den- 
sities of states, Eq. (Al). This is possible because the 



then given by 



fluctuations correspond to excitations of virtual pairs of 
particle-hole pairs. The particle-hole densities of states 
can be calculated analytically as functions of temperature 
and magnetization [see Appendix A] . We can write down 
the expression for the fluctuation contributions to the 
Ginzburg-Landau coefficients, a a and /3g in terms of in- 
tegrals over the modified particle-hole densities of states 
and their derivatives with respect to the magnetization 
[see Appendix B] . The fluctuation corrections an and /3fl 
are computed analytically at low temperatures and nu- 
merically over the full temperature range. In this work, 
the phase diagrams are calculated using only the mean- 
field contribution to the 6th order coefficient 7, since it is 
a higher order term in the Ginzburg-Landau expansion 
and since the fluctuation corrections to 7 are extremely 
difficult to compute. At low temperatures, the fluctua- 
tion contributions to the coefficients are given by 



a fl ~ -A(l + 21n2) 3 4 , 
A^l + ln-Vg 6 , 



(17) 



with A = [16v / 2]/[3(2tt) 6 ]. Note that here and in the 
following, g is given in dimensionless units. The ln(T//i) 
dependence of j3 is a remnant of the A/ 4 ln[Af 2 + (T/ fi)~] 
term in the free energy of Belitz et a/P^ 



B. Spiral 

Next, we calculate the coefficients of the Ginzburg- 
Landau expansion allowing for a spatial modulation of 
the magnetic order. We restrict our consideration to a 
single planar spiral. We exploit the fact that the free 
energy is a functional of the mean-field dispersion in the 
presence of spiral magnetic order. 



1. Mean-Field Dispersion in the Presence of Spiral 
Magnetic Order 

First, we determine the mean-field dispersion in the 
presence of spiral magnetic order. Let the directrix of 
the spiral wave-vector point along the z direction. For a 
planar spiral, the magnetization vector is then restricted 
to lie in the xy plane, M(r) = Af(cosQ • r,sinQ • r, 0). 
Note that the Hamiltonian (fTl) does not favor a particular 
direction of the spiral. The mean-field Hamiltonian is 



+ 9M 2 



where 



4>l 



k+Q/2,T' ^k-Q/2,4. 



(18) 



(19) 



Diagonalization of this Hamiltonian leads to the mean- 
field dispersion 



c k+i 




+ ( 5 M)2. (20) 



For the case of a quadratic dispersion, this reduces to 



(21) 



We see that the spiral wave-vector enters the energy dis- 
persion as an angle-dependent magnetization. Note that 
we have absorbed a Q 2 term into the chemical potential. 



2. Mean-Field Ginzburg-Landau Coefficients 

We will now make use of this mean-field electron dis- 
persion in the presence of a spiral in order to determine 
the Ginzburg-Landau expansion coefficients. We Taylor 
expand the free energy of a spiral in powers on magne- 
tization M and wave- vector Q, keeping the terms up to 
6th order, 

T[M,Q] = (a + /3iQ 2 + 7 iQ 4 )Af 2 

+ (/3 + 72<3 2 )M 4 + 7^ 6 , (22) 

where the coefficients ,$1,71,72 remain to be determined. 
A useful simplification at this stage is to rescale the spiral 
wavevector according to Q — > gQ, so that it has the 
same dimensions as M. In this way, (3 and /3\, and 7, 71 
and 72 have the same dimensions. Let us first consider 
The mean-field contribution is given by 



/?i,mf — 2—5 



k Q 

k F Q 



,( 3 ) 



Since T -C /J., derivatives of Fermi functions are strongly 
peaked around the Fermi energy and we can set |k| = hp 
in the scalar product which leads to a simple angular 
weight, 



k Q 



cos 9, 



where 9 is the angle between the vectors k and Q 
ter carrying out the angular integral we obtain ^i,mf 
|/?mf- Similarly, we obtain the proportionalities 7i,mf 
§7mf and 7 2 ,mf 



Af- 



7mf- 



(> 



3. Fluctuation Corrected Ginzburg- Landau Coefficients 



We introduce a d-wave spin-nematic order parameter 



Now, we proceed to analyze the fluctuation corrections 
to the expansion coefficients in Eq. ( 22 ) . As in the eval- 



uation of the mean-field coefficients it turns out that the 
fluctuation contributions to the expansion coefficients in 
Q are related to those of the uniform ferromagnet by an- 
gular factors which are identical to those found in the 
mean-field case. 

The fluctuation corrections to the free energy are given 
by an integral over momenta ki , . . . , k4 of a kernel that 
explicitly depends on each of the momenta through the 
mean- field dispersion (21 ). The fluctuation contributions 



to the Ginzburg-Landau coefficients are obtained by dif- 
ferentiating Eq. ^ with respect to M and Q. First we 
differentiate the kernel with respect to dispersion and 
then the dispersion with respect to M and Q. For exam- 
ple, the fluctuation contribution to the M 2 Q 2 coefficient 
is given by 



d 2 ^ 



3M 2 dQ 2 



Q=0,M=0 



We use two important facts in order to calculate this, 

i. that the free energy is a functional of dispersion and 

ii. that the spiral wave vector enters the mean-field dis- 
persion Eq. (21| like an angle dependent magnetization. 



The dispersion enters for each of the momenta k; in 
the momentum sum in Eq. ([8]). Differentiating with 
respect to Q 2 , therefore, will bring down factors of 
(kj • Q/(k F Q)) 2 , each of which will contribute with an 
angular factor as in the mean-field case. This leads to 
the proportionality /3i.h ps §/3fl- Combining this with the 
identical result for the mean-field contribution we obtain 
/?i fa |/3. When the proportionality between all of the 



coefficients is taken into account the free energy ( 22 ) be- 
comes 



T[M,Q] = ( a +-pQ 2 + -jQ*)M 2 
3 5 

+(f3 + 7Q 2 )M i +jM 6 . 



(23) 



N 



d * = p-(^-fc 2 )«sin 2 0cos(2<£), (24) 



where c?k is the d-wave distortion. The spin-nematic or- 
der parameter looks like a magnetization order param- 
eter weighted by an angular factor. It corresponds to 
Fermi-surface distortions which have opposite signs for 
spin-up and spin-down electrons. The net magnetiza- 
tion however, vanishes since the volumes of the distorted 
spin-up and spin-down Fermi surfaces are the same (see 
Fig. |ljc)). As we will see later, it is straightforward to 
generalize our final results to spin-nematic states with 
different symmetries. 

Turning on a fictitious field hjy conjugate to the ne- 
matic order parameter, the dispersion is modified to 
e£ = £k — a h pf sin 2 8 cos(2</>). Next we expand the gener- 



ating function <p in powers of Hn , 



tp[h 



N 



a*h 2 



N 



(25) 



where the superscript tp is used to distinguish coefficients 
of the generating function from those of the Ginzburg- 
Landau expansion. 



1. Mean-field contributions to the coefficients of the 
generating function 

The mean-field coefficients of the generating function 
are given by 



<f = 9 2 {dl)'£nW(e k ) = {d 2 k )(a MF -g) 



MF 



4! 



9 4 <4>5>( 3 )(e k ) = <4)/3 



MF 



- 5 6 (4)£n( 5 )(e k ) = (4)7, 



(26) 



C. Other Instabilities 

Other order parameters, not driven in mean field, are 
slightly more difficult to analyze. This is because we can- 
not decouple the interaction in those channels. Instead, 
we introduce a field conjugate to the order parameter 
and calculate the generating function. We show how the 
coefficients of the generating function are related to the 
corresponding Ginzburg-Landau coefficients of the uni- 
form ferromagnet. A Legendre transform of the gener- 
ating function recovers the free energy. We outline this 
procedure for the case of a spin nematic. 



where (...) = ^ J^dOsmO f^dcf)... denotes the an- 
gular average. The coefficients are proportional to the 
corresponding uniform Ginzburg-Landau coefficients; the 
constants of proportionality are angular averages of pow- 
ers of the nematic distortion. The resulting integrals 
are easy to calculate, yielding (g? 2 .) = ^, (d k ) = j^g, 
and (d^) = 5555 • Note that the term linear in g in the 
quadratic mean-field coefficient is absent since the inter- 
action is local in position space while the nematic distor- 
tion is local in momentum space, leading to a vanishing 
weight in the nematic channel. 
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2. Fluctuation Correction to Coefficients of the Generating 
Function 

As we found in the case of the spiral, the fluctua- 
tion corrections to the coefficients in the nematic ex- 
pansion are related to those of the uniform ferromagnet 
by the same angular averages as the mean-field coeffi- 
cients. For example, let us consider the fluctuation con- 
tribution to the ./V 4 coefficient. When differentiating Q 
with respect to four times, this brings down terms 
like (d^d^d^d^) , where |k«| ~ kp since derivatives 
of Fermi functions are peaked around the Fermi energy. 
Angular averages of this type are potentially more com- 
plicated as the directions of different k's are not inde- 
pendent. However, the fact that the dominant contri- 
bution comes from the particle-hole pairs with momenta 
|ki— k 3 | = |k 2 — k4| « 2kF leads to a tremendous simplifi- 
cation. Within this approximation, ki , k 2 , k 3 and k4 are 
either parallel or antiparallel to one another, rendering 
(d^d^d^d^) w (dfy- Similar arguments hold for other 
types of terms that appear in the expansion. Thus, to 
leading order, the same proportionality holds as for the 
mean-field coefficients and consequently, the generating 
function is given by 



<p[h N ] = {di)(a-g)h% + (di)0h% 



(27) 



3. Ginzburg- Landau expansion of Spin Nematic 

In order to obtain the Ginzburg-Landau expansion of 
the free energy of the d-wave spin nematic, we perform 
the Legendre transform to leading order, 



F[N] = <p[h N [N]] + h N N 
dip 
dh~ 



-N. 



N 



(28) 



Carrying out this transformation, we obtain the free- 
energy expansion in powers of the nematic order param- 
eter N, 

F[N] = -{dl)(a - g)N 2 + (di)f3N 4 + (d^N 6 . (29) 

We point out that the derivation of the free-energy func- 
tional for the spin nematic is general and not constrained 
to a particular symmetry of the distortion. To obtain 
the free energy of a p-wave spin nematicpiHH] we simply 
replace the d-wave distortion dk by the p-wave angular- 
weight pk « cos 9. This leads to slightly different angular- 
averages, (p k ) = |, (p k ) = |, and (p k ) = | and hence 
to slightly different coefficients in the Ginzburg-Landau 
expansion. 



D. Deviation from Free-Electron Dispersion 

Our approach also enables us to analyze the modifica- 
tions to the phase diagram in the presence of a dispersion 



that slightly deviates from the free-electron dispersion, 
£k = \ + <5ek- We plug this expression into the general 
mean-field dispersion in the presence of a spiral, Eq. (20 1, 
to obtain 



k 2 



1. 



-5e k -<7\ Q ■ (k + -V5e k 



(ghiy. (30) 



In order to calculate the corrections to the Ginzburg- 
Landau coefficients, we first differentiate the free energy 
with respect to the dispersion and then the dispersion 
with respect to M and Q. Finally, we expand the re- 
sulting Ginzburg-Landau coefficients in powers of Se^, 
assumed small. The free energy is now given by the sum 
of Tm,Q (23) and the additional contribution 
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where (...) denotes an angular average and we have as- 
sumed that the deviation <5ek from the free-electron dis- 
persion is such that the odd-power angular averages give 
zero. Mixing between coefficients at different total order 
in M and Q now occurs, since the angular distortion 
enters in both spin-symmetric and spin-antisymmetric 
ways, as opposed to the spin-antisymmetric Fermi- 
surface distortion of spiral or spin-nematic states in the 
isotropic case. 

Similarly, the free energy of the spin-nematic state is 
now given by the sum of T\N\ (29) and an additional 
term given by 



8F[N] 
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(32) 



IV. PHASE DIAGRAM 



In the previous section we have shown that the free 
energy is a functional of the mean-field dispersion. We 
have used this fact to develop the Ginzburg-Landau ex- 
pansions for the uniform ferromagnet and the spiral fer- 
romagnet. The derivation of the free-energy of the spin 
nematic is more complicated since the point interaction 
has no weight in the spin-nematic channel and a mean- 
field decoupling is not possible in that case. Instead, we 
added a field conjugate to the nematic order parameter, 
calculated the generating function and then performed a 
Legendre transform to obtain the Ginzburg-Landau ex- 
pansion in terms of the spin-nematic order parameter. 
Finally, we considered small deviations from the isotropic 



free-electron dispersion and derived the resulting correc- 
tions to the free-energy functionals. 

We have shown how a very useful simplification oc- 
curs in the low temperature regime; all spiral and spin- 
nematic coefficients can be related to those of a uniform 
ferromagnet, which we have calculated analytically at low 
temperatures. 

We now use the Ginzburg-Landau functions developed 
above to construct the phase diagram as a function of 
temperature, T, and renormalized interaction strength, 
g. We minimize the Ginzburg-Landau free energy with 
respect to the order parameter(s) and compare the free 
energies of different phases. We show how quantum fluc- 
tuations stabilize the spiral and spin-nematic phases, nei- 
ther of which are favored in mean-field theory. The ef- 
fect of a small anisotropic correction to the free-electron 
dispersion on the topology of the phase diagram is also 
investigated. 



A. Mean-Field Phase Diagram 

First, we consider the mean-field phase diagram. From 
Eq. ( 16 ) we see that in the low-temperature regime, the 
Pmf, 



quartic coefficient, /3mf, is positive. The second order 
phase transition between paramagnetic and uniform fer- 
romagnetic state happens when the quadratic coefficient, 
QMF] changes sign. The mean-field phase diagram is 
shown in Fig. [2j Neither spiral, nor spin-nematic states 
are favored in mean field. 



1. Uniform, Ferromagnet 

Before investigating how fluctuations may favor the 
formation of new phases, we first investigate their effect 
on the uniform ferromagnet. From Eq. (17) we see that 



the fluctuations provide a negative contribution to the 
Ginzburg-Landau coefficients. Ferromagnetism is thus 
favored for lower values of the interaction strength g than 
in the mean-field theory. This becomes evident if we con- 
sider the line of second order transitions a = 0. 
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FIG. 3: (Color online) Phase diagram of the uniform ferro- 
magnet, including quantum fluctuations. Below the tricritical 
point (shown as circle), quantum fluctuations drive the tran- 
sition first order. 
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FIG. 2: (Color online) Phase diagram of the uniform fer- 
romagnet in mean-field theory. The transition between the 
uniform ferromagnet and the paramagnet is always second 
order. 



B. Fluctuation-Corrected Phase Diagram 

Next, we include quantum fluctuations in our analysis 
and allow for the generation of new phases that were not 
present in the mean-field theory. 



In the presence of fluctuations, the quartic coefficient 
(5 inevitably becomes negative for low enough temper- 
atures due to the ~ ln(T//x) divergence. For /3 < 
the paramagnet-to-ferromagnet transition becomes first 
order. The line of first order transitions is given by 
/3 2 = Aaj (the condition for degeneracy of the minima 
of the Ginzburg-Landau function). The location of the 
tricritical point, at the intercept of the first-order and 
second-order lines, is found at T* = 0.24^. This is in 
good agreement with previous numerical calculations.^ 
The occurrence of first-order transitions at low tempera- 
tures has been observed in numerous experiments!*^"^ 



2. Fluctuation- Driven Spiral 



From Eq. (22 1 we see that the Q 2 M 2 term favors non- 
zero Q for Pi < 0. The particular relationship between 
coefficients that is found in the free-electron case implies 
that this occurs when (3 < 0, i.e. the spiral first forms at 
the tricritical point where the transition into a uniform 
magnet would have become first order. The phase dia- 
gram showing the instability towards the formation of a 
magnetic spiral is shown in Fig. [4] We now derive it from 
the Ginzburg-Landau functional. 

In the case of the free-electron dispersion we have 
shown that fix = 2/3/3. This results in the formation 
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of a spiral state below the tricritical point. Minimizing 
the free energy Eq. (23) with respect to Q, we obtain the 
optimal wave vector 



2 



5 



(33) 



After substituting this value of Q back in ( 23 ), we obtain 
the free energy at the optimal wavevector as a functional 
of M, 

T Q [M] = a Q M 2 + B Q M 4 + 7Q M 6 , (34) 
27 7 



a Q = 
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(a) Spiral-to- Paramagnet Transition. In principle, 
there are two ways in which the system can make a 
transition from the paramagnet into a spiral state: 

i. A second order transition in M, along which 
M = 0. This line is given by aQ = 0, and the necessary 
condition for its existence is that Bq > 0. 

ii. A first order transition in M, along which M 
jumps from zero to a finite value. This transition 
happens along the line Bq = Accq^q, as long as Bq < 
and an > 0. 




FIG. 4: (Color online) Phase diagram of the itinerant fer- 
romagnet allowing for the possibility of spatially modulated 
ferromagnetism. Quantum fluctuations drive the formation 
of a spiral phase which sets in below the tricritical point and 
pre-empts the first-order transition between the uniform fer- 
romagnet and the paramagnet. 

Since we have already established that we can have a 
spiral state only for Bi < and hence (3,Bq < (fol- 
lowing from the proportionalities of coefficients) , we rule 



out the first possibility and conclude that the transition 
from the paramagnet into the spiral ferromagnet must be 
first order in M (and also first order in Q, according to 
Eq. (33)). Substituting 0£q, Bq and 7q from Eq. (34), 
the equation for this line becomes aj = ||/3 2 . This tran- 
sition pre-empts the transition from the paramagnet into 
the uniform ferromagnetic state [see Fig. [4] . 

(b) Uniform Ferromagnet-to- Spiral Transition. Next, 
we wish to determine the boundary between the spiral 
phase and the uniform ferromagnet. In principle, this 
may occur either discontinuously or smoothly. In the 
case of the free-electron dispersion, it occurs via a Lifshitz 
transition where Q goes continuously to zero. The value 
of magnetization that minimizes the free energy Tq is 
given from Eq. (34) by 




c*7 



5 

27 



(35) 



Substituting this into Eq. ( 33 1 for Q, we find that the Lif- 
shitz line coincides with line a — 0. The magnetization 
M varies continuously along this line. 
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FIG. 5: (Color online) The evolution of the order parameters 
M and Q in the spiral phase. As we approach the tricritical 
point, the jumps in M and Q along the line of first order 
transitions become smaller. At the tricritical point M = Q = 
0. At the Lifshitz transition between the uniform ferromagnet 
and the spiral ferromagnet, Q goes continuously to zero while 
M remains finite and behaves smoothly. 

The evolution of the order parameters M and Q in the 
vicinity of the first order transition from the paramagnet 
into the spiral state are plotted in Fig. [3J We see that 
the jumps in M and Q get smaller as we approach the 
tricritical point. 

Previous analyses^ determined the phase diagram of 
spiral and uniform ferromagnets (they did not consider 
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a spin nematic phase) by brute force numerical (Monte- 
Carlo) evaluation of the fluctuation corrections given by 
Eq. ^ for given g, T and M,Q. We, instead, evaluate 
the phase diagram within a Ginzburg-Landau expansion 
and obtain an analytical approximation at low tempera- 
tures. The agreement between the two methods is good 
in the vicinity of the tricritical point where the expansion 
is controlled. In addition, we were also able to determine 
the location of the boundary between the uniform and 
spiral ferromagnets as well as the nature of this transi- 
tion. 



3. Fluctuation-Driven Spin Nematic 

Finally, we determine the region of the phase diagram 
where the d-wave spin-nematic phase has the lowest free 
energy. By analyzing the free energy Eq. ( 29 ) , we find 
that for temperatures below T = 0.02/i, 



the transition 

from the paramagnet into the spiral state is pre-empted 
by a transition into a spin-nematic state. The first or- 
der transition line between the paramagnet and the spin 
nematic is given by the equation 



ah- 



(36) 



From the evaluation of this equation for spin-nematic 
states with d- and p-wave symmetry we also find that the 
instability to the formation of the d-wave spin nematic 
occurs at slightly higher temperatures and is therefore 
favored. This however might change with dimensional- 
ity, the form of the electron dispersion, or the range of 
the interactions. 
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FIG. 6: (Color online) Phase diagram of the itinerant ferro- 
magnet, allowing for the possibility for the formation of spiral 
and spin-nematic phases. At temperatures which are about 
an order of magnitude smaller than the temperature of the 
tricritical point, a d-wave spin nematic forms between the 
spiral ferromagnet and the paramagnet. 

Comparing the free energies of the spiral and spin- 
nematic phases, we find that the spin-nematic state pen- 



etrates into the region where the spiral state was pre- 
viously favored. The details of the transition between 
the spiral and spin-nematic phases are potentially very 
interesting but hard to analyze. Introducing phase slips 
into the spiral restores translational symmetry and ren- 
ders the phase nematic. Whether this is indeed the spin- 
nematic phase identified here, or something more exotic^ 
is not clear. 

The phase diagram, including the spin-nematic state is 
shown in Fig. [6j Note that this phase diagram is plotted 
on a logarithmic scale, since the spin-nematic state onsets 
at a temperature which is an order of magnitude lower 
than the temperature of the tricritical point where the 
spiral forms. 

In summary, quantum fluctuations have generated a 
coupling in the spin-nematic channel and stabilized a 
spin-nematic phase. This is similar to the mechanism 
by which a super condu cting state is stabilized in spin- 
fluctuation theoryj2Sl23 anc } we anticipate that our ap- 
proach can be applied to study superconductivity as well. 
We emphasize that the quantum order-by-disorder ap- 
proach incorporates charge fluctuations on the same foot- 
ing as spin fluctuations (see Section II C) . As pointed out 
by Chubukov and MaslovJ^ charge fluctuations are es- 
sential to mediate the formation of a spin-nematic state. 



C. Phase Diagram for Anisotropic Dispersion 

Changing from a free-electron dispersion to a band dis- 
persion can have a profound effect upon the magnetic 
phase diagram. In the extreme, it can lead to nest- 
ing and the formation of modulated (anti-ferromagnetic) 
phases even at mean field. Here, we consider the effect 

,2 

of a weakly anisotropic dispersion £k = ^ + ^k; with 
<5ek = 8 cos (40) sin 9. This deformation makes the dis- 
persion more tight-binding like. By changing the subtle 
balance between kinetic energy and fluctuation correc- 
tions the topology of the phase diagram, the regions oc- 
cupied by the different phases, and the nature {e.g. first 
or second order) of the transitions are altered. 



1. Uniform Ferromagnet 

First, we investigate the changes to the phase diagram 
of the uniform ferromagnet. Summing Eqs. ( 15 ) and ( 31 ) 
evaluated at Q = 0, we arrive at the following expression 
for the free energy: 



T = c\M 2 + f3M 4 + 7 M e 
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As in the case of the isotropic k 2 dispersion, we find 
that the transition between the uniform ferromagnet and 
the paramagnet is continuous at high temperatures and 
becomes first order at low temperatures due to a sign 
change of $. The line of second order transitions between 
the uniform ferromagnet and the paramagnet is given by 
a = while the line of first order transitions is given by 
j3 2 = 4S7. The effect of the anisotropic correction to the 
dispersion is to slightly shift the locations of the phase 
boundaries, e.g. the temperature of the tricritical point 
is reduced to T* = 0.225/; (see Fig. [7}. 

As we will see in the following, the effects of the 
anisotropy on the formation of the spiral and the na- 
ture of the transitions to the spiral ferromagnet are more 
interesting. 



2. Fluctuation- Driven Spiral 

Our analysis of the spiral phase follows the same steps 
as in the case of the free-electron dispersion in Section 
IV. 2. The resulting expressions are lengthy and not par- 
ticularly revealing in themselves. Therefore, we simply 
outline the main steps. The free energy of the spiral state 
is the sum of Eqs. (23) and (31 1, and is given by 



F[M, Q] = (a + /?i(Q)Q 2 + 71 Q 4 )M 



•721 



M + 7M , 



with a, (3, and 7 defined in Eq. (|37|), 71 = 71 

72 = 72 = li and 



(38) 



? 7, 
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This free energy now depends upon the direction Q = 
Q/Q and is no longer invariant under rotations of the 
spiral. This is the consequence of the anisotropic dis- 
persion which breaks the continuous rotation symme- 
try It turns out that for the particular anisotropy 
ckk = S cos (4<^>) sin 6, the free energy is minimized for 
spirals with Q along the z-axis. 

Notice that the proportionality between coefficients 
found in the case of the free-electron dispersion is bro- 
ken by the anisotropic dispersion. For example, the co- 
efficient of the term Q 2 M 2 is no longer proportional to 
that of the M 4 term. This broken symmetry changes the 
topology of the phase diagram and changes the nature 
of the transition between the spiral ferromagnet and the 
paramagnet (see Fig. 

(a) Optimum wavevector. The optimum wavevector 
is obtained by minimizing the free energy Eq. ( 38 1 with 



respect to Q for a given magnetization M. In this way, 
we obtain Q = Q[M]. Substituting this value of Q back 
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FIG. 7: (Color online) Modification to the phase diagram of 
the spiral state with a weakly anisotropic dispersion. The on- 
set of the spiral no longer coincides with the uniform tricrital 
point (circle). Instead, the spiral forms at a slightly higher 
temperature (square) and pre-empts a portion of the continu- 
ous transition between the uniform ferromagnet and the para- 
magnet (thin dashed line) as well as the first-order transition 
(thick dashed line). Note that because of the anisotropy also 
the nature of the spiral-to-paramagnet transition changes. At 
higher temperatures M now behaves continuously while at 
low temperatures the transition is first order in M as in the 
isotropic case. 



into Eq. (38) we obtain the free energy at the optimum 
wavevector as a function of M, 



PqW] = &qM 2 + p q m 4 



7q 



(40) 



where olq, /3q and are the appropriate generalizations 
of the free-electron forms given in Eq. ( 34 ) . 

(b) Paramagnet-to-Spiral Transition. As discussed in 
the case of the isotropic free-electron dispersion, in prin- 
ciple, there are two ways in which the system can make 
a transition form the paramagnet into the spiral state: i. 
a second order transition in M (which occurs along the 
line olq = for /3q > 0), and ii. a first order transition in 
M (which occurs along the line /3q = AcyqJq for /3q < 
and otQ > 0). For the free-electron dispersion only the 
latter possibility occurred. Now that we have broken the 
proportionality between Ginzburg-Landau coefficients by 
allowing for an anisotropic dispersion, both of the possi- 
bilities exist. 

i. The transition begins as second order in M (and first 
order in Q) (thin solid line in Fig. [7J extending below 
the square symbol). This line was not present in the 
phase diagram with free-electron dispersion, because the 
quartic coefficient j3 was negative in the region where the 
spiral existed. This line pre-empts a portion of the line 
of second order phase transitions between the paramag- 
net and the uniform ferromagnet (thin dashed line) and 
the line of first order transitions from the paramagnetic 
to the uniform ferromagnetic state (thick dashed line). 
The formation of the spiral phase is, therefore, no longer 
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coincident with the putative tricritical point of the uni- 
form ferromagnet (shown as circle in Fig. [7]) but sets in 
at higher temperatures as indicated by a square symbol 
in Fig. 

ii. The second order transition between the paramagnetic 
and spiral phases gives way to a first order transition at 
lower temperatures - shown as thick solid line in Fig. [7j 
(c) Uniform Ferromagnet-to- Spiral Transition. The 
boundary between the spiral and uniform ferromagnetic 
phases remains a Lifshitz transition, where the optimal 
wave vector falls continuously to zero and M is continu- 
ous. In the case of the free-electron dispersion, this has 
turned out to be coincident with the a = line. The sit- 
uation is not so simple when we allow for an anisotropic 
dispersion. While the magnetization remains continu- 
ous, within our numerical resolution, we cannot exclude 
that the derivative of M becomes discontinuous. The 
anisotropic dispersion has had two key effects upon the 
phase diagram. Firstly, the region occupied by the spi- 
ral phase has increased, and secondly, the onset of the 
spiral decoupled from the tricritical point of the uniform 
ferromagnet. 



3. Fluctuation-Driven Spin Nematic 

The free energy of a d-wave spin-nematic state in the 
presence of the distortion is given by the sum of Eqs. ( 29 ) 
and ( 32 1 . We can rewrite this in the form 



F[N] = & N N 2 + fi N N 4 + j N N 6 



(41) 
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The transition from the paramagnet into the d-wave spin 
nematic is first order and occurs when j3fj = ^czn^n- 
The boundary between spiral and nematic states is ob- 
tained by comparison of their free energies - although, as 
we stated above, the details of how this transition occurs 
may be subtle. When the dispersion is anisotropic, the 
spin nematic state occurs at higher temperatures than 
in the isotropic case. The spin nematic is stabilized 
since it becomes easier to redistribute the kinetic energy 
cost of forming a spin nematic when the dispersion is 
anisotropic. 

In summary, we have used the Ginzburg-Landau ex- 
pansion of the free energy, including the quantum fluctu- 
ations, to determine the phase diagram of the itinerant 
ferromagnet in the vicinity of a quantum critical point. 
First, we investigated the effect of quantum fluctuations 
upon the phase diagram of the uniform ferromagnet. Be- 
low a certain temperature (the tricritical temperature) 
the paramagnet-to-ferromagnet transition becomes first 



order. Next, we allowed for the possibility of a spa- 
tially modulated ferromagnetic phase. For temperatures 
lower than the tricritical temperature, it becomes ener- 
getically favorable to form a spiral state in between the 
paramagnetic and uniform ferromagnetic states. The 
putative first order transition between the paramagnet 
and the uniform ferromagnet is pre-empted by a transi- 
tion into this spatially modulated state. Further, we in- 
cluded the possibility of forming a d-wave spin-nematic 
state. This state is stabilized for temperatures, below 
T R3 0.1T* w 0.02/j, in between the paramagnetic and 
spiral phases. Finally, we have shown that in the pres- 
ence of an anisotropic dispersion, the topology of the 
phase diagram changes such that both spiral and spin- 
nematic phases occupy a larger region of the phase di- 
agram. Moreover, the onset of spiral order no longer 
coincides with the uniform tricritical point. 



V. CONCLUSIONS AND OUTLOOK 

We have shown how quantum fluctuations can lead to 
the formation of new phases in the vicinity of itinerant 
ferromagnetic quantum critical points. Quantum order- 
by-disorder not only provides an intuitive physical pic- 
ture for this process but identifies a general principle 7 
behind the phase reconstruction near quantum-critical 
points. 

Quantum order-by-disorder is familiar in condensed- 
matter^^ as well as in high-energy physics where it 
is referred to as the Coleman- Weinberg mechanism!^ In 
these familiar realizations, new ground states are stabi- 
lized by quantum fluctuations of a bosonic order parame- 
ter. What makes our approach new is the immediate con- 
nection between Fermi-surface deformations (associated 
with the onset of competing order) and the enhancement 
of phase space available for low-energy quantum fluctua- 
tions. 

Recently, it has been arguecP^ that the avoidance of 
naked quantum-critical points due to the quantum order- 
by-disorder mechanism can be understood within the 
AdS/CFT correspondence^^ between conformal field 
theories, describing critical condensed-matter systems, 
and gravity in Anti-de-Sitter space. In the gravity con- 
text, the quantum-critical state at finite temperatures 
corresponds with a Reissner-Nordstrom black hole in AdS 
space. It has been realized that such a black hole can 
become unstable at low temperatures and tends to col- 
lapse to a state with lower entropyP^ This entropic mo- 
tive underlies the quantum order-by-disorder mechanism. 
Experimentally, the measurement of entropic landscapes 
has proven a revealing probe of the phase reconstruction 
near quantum critical pointsPS 

The fermonic quantum order-by-disorder approach 
presented in this paper not only establishes the connec- 
tion to deformations of the Fermi surface, which are ac- 
cessible by various experimental probes, but also leads 
to relatively simple analytical calculations, based on self- 
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consistent second order perturbation theory. As such, 
it is more accessible than technically involved diagram- 
matic techniques ESHH1 The two approaches are formally 
equivalent; expanding self-consistently about a saddle 
point with the already established order re-sums selected 
series of diagrams that give rise to non-analytic correc- 
tions to the free energy. 

We have used this approach to investigate the 
fluctuation-driven phase reconstruction in the vicinity 
of the itinerant ferromagnetic quantum critical point in 
three spatial dimensions. This quantum-critical point 
is unstable towards the formation of spiral and spin- 
nematic states. Quantum fluctuations would render 
the transition between the uniform ferromagnet and the 
paramagnet first order. This first-order transition is pre- 
empted by a modulated, spiral ferromagnetic phase. At 
even lower temperatures, a d-wave spin-nematic state 
forms which is slightly favored over a spin-nematic with 
p-wave symmetry. It is sandwiched between the para- 
magnetic and spiral phases. In order to describe more 
generic experimental systems, we determined the topol- 
ogy of the phase diagram in the presence of an anisotropic 
electron dispersion. The regions of phase space occupied 
by both the spiral and spin nematic phases are enlarged. 
Moreover, the onset of spiral order and the putative tri- 
critical point of the uniform ferromagnet become decou- 
pled and the order of transitions is modified. 

Similar spin-nematic instabilities near to the itinerant 
ferromagnetic quantum critical point have been studied 
recently by Chubukov and Maslov 23 within an extension 
to Hertz-Millis theory in a spin-fermion model. The au- 
thors point out that the inclusion of the effects of charge 
fluctuations (through Aslamov-Larkin corrections) in ad- 
dition to spin fluctuations is crucial to mediate the forma- 
tion of a spin nematic. Our quantum order-by-disorder 
approach incorporates charge fluctuations on the same 
footing as spin fluctuations. The results of Ref. [23] are 
very much in accord with those that we present here; a 
spin nematic instability occurring out of the paramag- 
netic state is found to pre-empt the spiral phase. We 
point that the contact interaction used here as opposed 
to the finite-range interaction in Ref. [23] lowers the free 
energy of the spiral relative to the spin nematic. 

The quantum order-by-disorder approach can be ap- 
plied to a variety of systems and phases. Adding a small 
spin-orbit coupling to the Stoner model of magnetism 
had previously enabled us to explain the partially ordered 
phase of MnSPin terms of quantum order-by-disorder 
Work on the superconducting instability is in progress. 
It appears that our approach recovers the results of spin- 
fluctuation theory, thus revealing the link between spin- 
fluctuation theory used to treat superconducting insta- 
bilities, and the extensions of Hertz-Millis theory that 
explain the instability of quantum critical points to other 
types of order. 

There are several natural directions for developing our 
approach. The inclusion of superconducting instabilities 
is a priority - the nature of the superconducting phase 



where it overlaps with the spatially modulated magnetic 
phases raises the intriguing possibility of spontaneous, 
fluctuation- driven, modulated superconductivity. The 
calculations described in this paper are in three dimen- 
sions. Extending them to two dimensions, where the non- 
analytic effects of fluctuations are even stronger, is an 
important step. Finally, the role of fluctuations near to 
anti-ferromagnetic quantum critical points may also be 
susceptible to analysis using our methods. 
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Appendix A: Free energy of the uniform 
ferromagnet in terms of modified particle-hole 
densities of states 

The fluctuation corrections to the free energy are given 
by a high dimensional integral over momenta k 1; . . . , k4 
and correspond to excitations of virtual pairs of particle- 
hole pairs of opposite spin and equal and opposite mo- 
menta. It is therefore possible to rewrite the regularized 
fluctuation corrections J-^ Q as a lower dimensional in- 
tegral over modified particle-hole densities of states, 



CT=±1 Jq.ei, 



Ap CT (q,e 1 )p- J (-q, eg) 
ei + (2 



(Al) 



dr. 



where we have defined J q := J an d f e '■= f° 
The modified particle-hole densities of states as a func- 
tion of momentum q and energy e are given by 



p°(q,e) = / n(e£_ f )<5(e-e£ +f +e£_ f ), (A2) 
Ap a (q,e) = /n(e£_a)n(eJ +t )*(e-eJ +? +<_ f ), 

Jk 

and are related to the particle-hole density of states as 
pph = Ap CT — p a . This form of the fluctuation correc- 
tion will prove useful in our subsequent evaluation of the 
Ginzburg-Landau expansion of the uniform ferromagnet, 
since the modified particle-hole density of states of the 
uniform ferromagnet can be calculated analytically. This 
leads to a tremendous simplification of the fluctuation in- 
tegral. 

The modified particle-hole densities of states are func- 
tions of the magnetization M , which enters through the 
dispersion e£ = ek — agM of the uniform ferromagnet. 
We want to Taylor expand Eq. ( Al ) with respect to M. 



In doing so we will require the derivatives of p a and Ap a 
with respect to M. However, since in p" and Ap a the dis- 
persion only enters for either spin up or spin down (and 
not both) we can relate the derivatives with respect to 
M to derivatives with respect to the chemical potential 



14 



/<■, 

d\ 4 Ap°{<l^)\ M=0 = (<jgydlAp°( qi e)\ M=0 

= (ag)%Ap(q,e). (A3) 

Now let us derive explicit expressions for Ap = Ap a \ AI=0 
and p — p a \ M=0 and their derivatives. 

1. Evaluation of p(q,e) and Ap(q, e) 

For an isotropic free-energy dispersion, = |fc 2 the 
angular integrals in Ap(q, e) and p(q, e) are easy to com- 
pute in three dimensions since the only angular depen- 
dencies enter through the volume element and the scalar 
product k • q = kq cos 8. The remaining one-dimensional 
integrals over k are elementary. The resulting modified 
particle-hole densities of states are given by 




0.05 0.1 0.15 0.2 0.25 



Tjp 

FIG. 8: (Color online) Fluctuation contributions to a: com- 
parison of numerics (solid line) with leading low-temperature 
analytical dependence (dashed line). 



Ap(<?,e) 
where, 



1 T 



(2tt)2 q 
1 T 
(2^9 



m M +e -tV , w;-m) j (A4) 
1 — eT \ J 



l + e"T V / 



(A5) 



For the derivatives of the modified densities of states 
with respect to the chemical potential wc obtain 



<9>(<7,e) 
dlAp(q,e) 
where n(x) 



1 1 



d^n [<t>+(e, q)], (A6) 
d$-»{n[r(*,q)] n[^(e,q)}}, 



(27r) 2 q » 

1 1 
(2tt) 2 g"" 

l/[eS x ~^/ T +l] denotes the Fermi function. 



Appendix B: Fluctuation contributions to a and j3 

Here, we outline the evaluation of the fluctuation con- 
tributions to the Ginzburg-Landau coefficients a and j3 
of the uniform ferromagnet, 



1 a 2 j- fl 



2 dM 2 



and /3fl 



1 a 4 Jfl 



M=0 



4! <9M 4 



(Bl) 



M=0 



The modified particle-hole densities of states and their 
derivatives with respect to the chemical potential were 
calculated in Appendix [A] For brevity, we define the in- 
tegrals 



9l,Ap(q, ex)dlp(q, e 2 ) 



q, e i. e 2 



ei + £2 



(B2) 



The fluctuation contribution to the free energy at zero 
magnetization = Fa\ M=0 is given by 



(o) 



l T 

■JQ.0- 



(B3) 



We evaluated the integral Jq,o numerically for finite tem- 
peratures and analytically at T = 0. The analytical 
calculation at T — correctly reproduces the result of 
Abrikosov and Khalatnikov.321 By Taylor expanding the 
fluctuation corrections to the free energy, we obtain the 
fluctuation contributions to the Ginzburg-Landau coeffi- 
cients, 



a a = 2 5 4 ]T(-1) 



i=0 



Ji.1-i 



/3fl 



tB-WJW*. ( B4 ) 

i=0 v 7 



Some of the integrals Jj j are difficult to evaluate numer- 
ically, since at very low temperatures the higher deriva- 
tives of the Fermi functions (which enter through the 
derivatives of modified particle-hole densities of states) 
are strongly peaked around the Fermi energy and rapidly 
change sign. 



Evaluation of a 



From Eq. ( B4 ) we see that we need to calculate three 



terms Jo, 2,^2,0 and Ji.i- In principle, J0.2, J2,o are more 
difficult to calculate numerically (but possible) since they 
contain first derivatives of Fermi functions which are 
strongly peaked. However, we can reduce the number 
of terms that we need to calculate by collecting some 



together to re-express them as derivatives of J^ 
respect to p. For example, 



(0) 



with 



(o) 



V(J ,2 + 2Ji,i+ J 2 



(B5) 
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Using these relations we can rewrite the fluctuation con- 
tribution to a as 



a a = 2<? 4 (J , 2 - 2J M + J 2>0 ) 

„2 



_ 3 fl 2r(0) 

~ 2! M fl 



(B6) 



The remaining integral Ji x is easy to evaluate numer- 
ically. The temperature dependence of aa is shown in 
Fig.|8] At zero temperature we obtain 



a fl ~-A(l + 21n2) ff 4 , 



where A 



3(277) 6 ■ 



2. Evaluation of /3 



(B7) 



Fermi functions they are even more difficult to evaluate 
numerically. We collect some of the terms together by 
noting that 



to get 



a _ 3 ^4 t-(O) 



"(Jl,3 + ^3,l)- 



(B8) 



(B9) 



In this way we avoid the explicit calculation of integrals 
Jo, 4 and Ja,o, which contain third derivatives of Fermi 
functions. We can further simplify by noting that 




FIG. 9: (Color online) Fluctuation contributions to j3: com- 
parison of numerics (solid line) with leading low-temperature 
analytical dependence (dashed line). 

In order to evaluate /3a we need to evaluate five terms 
in Eq. (B4). As they contain higher order derivatives of 



dl Ji,i = J 3 ,i + 2J 2 , 2 + Ji- 



(BIO) 



If we re-express 
we obtain 



Ji >3 in Eq. (B9| using Eq. (BIO I, 



/8fl 



(Bll) 



We have already calculated the functions J 7 ^ and Ji,i, 
when we evaluated aa. Both are smooth functions and 
the numerical evaluation of the derivatives with respect 
to /i is trivial. The leading temperature dependence of 
/3fl comes from the J 2 ,2 integral, which diverges as T —> 0. 
In this limit 



ft, - A^l + ln-]//. 



(B12) 



The logarithmic divergence with temperature arises from 
particle-hole pairs with momenta q ~ 2kp. The good 
agreement between our numerical and analytical results 
in the low-temperature regime is shown in Fig. [9j 
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